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Drumhead  Problem 


Given  a  sample  of  a  drums  sound,  attempt 
to  calculate  the  amplitudes  of  the  modes 

-time-series  at  location  of  microphone 

Related  to  a  famous  problem  posed  by 
mathematician  Mark  Kac  (1966)  asking: 
“Can  One  Hear  the  Shape  of  a  Drum?” 

Lead  to  idea  of  iso-spectral  drums 


Isospectral  Drums 

Drums  with  differing  boundaries  that  have 
identical  k_n  values  -  so  they  sound  alike! 
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Why  Measure  the  Amplitudes? 

Being  able  to  measure  modes  strengths 
suggests  dynamics  about  the  system. 

By  measuring  “windows”  in  time  that 
overlap,  the  time-dependence  of  the 
amplitudes  can  be  seen. 

Energy  conservation  -  once  a  mode  is 
excited,  where  does  the  energy  go? 

Lead  to  prediction  of  amplitudes  beyond 
the  time  window  (forecasting). 


Outline 
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Methodology 

Three  key  issues: 

-  Delta  Function 

-  Solution  Architecture 

-  Degeneracy  of  States 

Midn  1/C  Grant  Hundley 

-  Developed  a  drum  simulator 

-Working  on  numerical  scheme  to  extract 
amplitudes  from  simulated  time-series 


Delta  Function 
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Delta  Function  Assumptions 
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Wavelength  is  set  from  the 
window  and  is  symmetric 
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Define:  Delta  and  Epsilon 
Should  these  conditions 
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Delta  Function 
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assuming  m  is  free  (real) 

assuming  n  is  even  amd  m  is  an  integer 
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Delta  Function  ~  Sinc(pi*(n-m)) 


The  Anm  function  shown  below  shows  its  clear  approximation  to  6nm  when  (n,  m)  are  integers.  Also  shown  are  the  twin 
responses  at  (— m,  +m)  for  the  n  =  10  case. 


Delta  &  Epsilon  Function 


Delta  function  is  a  measure  of  the  “mixing 
between  the  eigenmodes. 

Slightly  different  Delta  fucntion  for  the 
cos*cos  term  (more  on  this  later). 

Epsilon  will  be  identically  zero  for 
symmetric  domains. 


Epsilon  Timing  Resolution  Function 


£(©m=10)  =  cos(jc*(1  0-x))/(jc*(1  0-x))  -  cos(ji*(1  0+x))/(ji*(1  0+x) 


Solution  Architecture: 
Solutions  to  space-time  problems 
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Nature  of  Eigenmodes  with  both 
space  and  time  in  solution 

Two  types: 

-  Coupled  -  for  each  k  eigenvalue  in  space 
there  exists  a  unique  co  eigenvalue  in  time 

•  Dispersion  relationship,  co(k),  for  most  differential 
equations  in  (x,t). 

•  Dispersion  relationship  is  generally  monotonic  and 
increasing. 

-  Decoupled  -  for  systems  not  well  motivated 
physically,  yet  can  be  described  with  a  space- 
time  basis  set  (river  bank  problem). 


Solution  Architecture 


At  one  location,  xl ,  sample  the  data  in  a  time-series. 
Project  out  the  sin  and  cosines. 
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Solution  Architecture 


•  Mixing  nature  of  Delta, Epsilon  matrix  can 
clearly  be  seen. 

•  No  nodes  can  exist  in  the  f(k*x)  matrix. 

•  Solve  for  Amplitudes  thru  matrix  inversion. 


Solution  Architecture 


Short-hand  notation: 
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Drum  -  Specific  Solutions 

Drum  problem:  membrane  stretched  over 
a  circular  boundary  (Dirichlet  be). 

Strike  the  drum. 

Use  a  microphone  to  record  the  time- 
series. 

Fourier  analyze  the  time-series  to  obtain 
amplitudes  for  each  com :  An 
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•  Spatial  modes  are  a  combination 

of  Bessel  function,  J(knr)  times  sin(n0) 
or  cos(n0) 

•  Temporal  modes  use 

sin(oont)  or  cos(cont) 

•  For  each  Bessel  function,  there  exists 
multiple  zero  crossings,  n’ 

•  kn  values  are  non-integerlike,  soconfail 
conditions  for  orthonormality 


Table  1:  Bessel  Function  Zero  Crossings 
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Table  2:  Bessel  Function  Zero  Crossings 
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Table  3:  Bessel  Function  Zero  Crossings 
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Degeneracy  of  States 
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Degeneracy  and  Sampling 

•  Due  to  symmetry  in  the  solution, 
degenerate  states  are  produced. 

•  Due  to  degeneracy,  there  is  a  2-1  ratio  of 
unknowns-knowns. 

•  Solution:  add  another  sample  location 
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Calculating  the  Amplitudes 
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Drum  Solution  in  Short-hand 


•  Number  of  sample  locations  “squares-off” 
the  degeneracy  matrix,  allowing  the 
system  to  be  solvable. 


Drum  Problems 

n=0  sin(n0)  term  needs  to  be  removed, 
breaking  the  2-1  ratio  to  less  than  2-1 . 

Test  method  against  drum  simulation,  with 
known  inputs  to  the  amplitudes,  An. 

Further  degeneracies  exist  due  to 
closeness  of  k-eigenvalues. 


Chesapeake  Bay  Problems 

No  known  analytic  solution  to  the  Bay 
(ie.  No  analytic  hints  as  to  any  degeneracy) 

At  a  given  samples  location,  calculate  the 
projections,  (am,pm)  for  a  range  of  (com,  T). 

Observe  the  patterns  of  oom  and  compare 
the  sequence  to  kn’s. 

Guess  the  dispersion  relationship  (map 
from  kn  to  com 


Degeneracy  of  States 

From  observables  at  frequencies  com, 
observe  projection  changes  as  the 
period, T,  is  changed. 

Select  an  appropriate  period,  T. 

Construct  Degeneracy  matrix  based  on 
best  guess  of  dispersion  relation  as  well 
as  k-eigenvalues  density. 


Timing  Resolution  Nature  of  Delta  Matrix 


ov-frequency 


Delta(oom,T) 


Delta((om,T) 
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Future  Plans: 


Run  thru  the  toy  model  (Drum). 

Add  source  terms  to  Chesapeake 
Bay  model. 

Add  Delta(spatial)  matrix. 
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Chesapeake  Bay  Problem 

Take  data  at  stations  around  the  Bay, 
collecting  time-series  of  vector  flows. 

How  many  stations  are  needed  to  provide 
enough  data  to  fully  calculate  the  modes? 


Image  Processing  of  the 
Chesapeake 
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Neumann  Mode  100  (350X1185) 
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Advantages 


Eigenmodes  fill  domain  (space),  suggest 
„„  r  future  behavior 

Monterey  Bay:  Surface  Currents 


